Inclusion of the first-order vector- and tensor-modes in the 
second-order gauge-invariant cosmological perturbation theory 
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Abstract 

Gauge-invariant treatments of the second-order cosmological perturbation in a four 
dimensional homogeneous isotropic universe are formulated without any gauge fixing. 
We have derived the Einstein equations in the case of the single perfect fluid without 
ignoring any modes. These equations imply that any types of mode-coupling arise 
' due to the second-order effects of the Einstein equations. 

(N : 

The second-order general relativistic cosmological perturbation theory has very wide physical moti- 
' vation. In particular, the first order approximation of our universe from a homogeneous isotropic one is 

revealed by the recent observations of Cosmic Microwave Background (CMB) by Wilkinson Microwave 
Anisotropy Probe [1], which suggests that the fluctuations of our universe are adiabatic and Gaussian at 
least in the first order approximation. One of the next theoretical researches is to clarify the accuracy of 
this result through the non-Gaussianity, or non-adiabaticity, and so on. To carry out this, it is necessary 
to discuss the second-order cosmological perturbations. 
' However, general relativistic perturbation theory requires delicate treatments of "gauges" and this 

situation becomes clearer by the general arguments of perturbation theories. Therefore, it is worthwhile 
to formulate the higher-order gauge-invariant perturbation theory from general point of view. According 
to this motivation, we proposed the general framework of the second-order gauge-invariant perturbation 
theory on a generic background spacetime; 2]. This general framework was applied to cosmological per- 
turbation theory[3j and all components of the second-order perturbation of the Einstein equation were 
derived in gauge invariant manner. The derived second-order Einstein equations are quite similar to the 
equations for the first-order one but there are source terms which consist of the quadratic terms of the 
linear-order perturbations. 

In this article, we show the extension of the formulation in Refs. [3J to include the first-order vector- 
and tensor-modes in the source terms of the second-order Einstein equation, which were ignored in 
Refs. [3j. 

As emphasized in Refs. [H [3J, in any perturbation theory, we always treat two spacetime manifolds. 
One is a physical spacetime A4\ and the other is the background spacetime A4q. In this article, the 
background spacetime M.q is the Friedmann-Robertson- Walker universe filled with a perfect fluid whose 
metric is given by 

g a b = a 2 (i]) (-{drj) a {drj)b + 7y \dx l ) a (dx J ) b ) , (1) 

where 7^ is the metric on maximally symmetric three space. The physical variable Q on the physical 
spacetime is pulled back to xQ on the background spacetime by an appropriate gauge choice X which is 
an point-identification map from Mq to M\. The gauge transformation rules for the pulled-back variable 
xQ, which is expanded as xQx = Qa + ^xQ + l^xQi are gi ven by 



Q - £( m Qo, ( yQ - { xQ - 2£ ^i) { xQ + { £ i(2) + ■£?(!)} Qo, (2) 



where X and y represent two different gauge choices, and £" 2 j are generators of the first- and the 
second-order gauge transformations, respectively. The metric g a b on the physical spacetime M.\ is also 
expanded as g a b — g a b + Aft. a fc + 4j-Z a fc under a gauge choice. Inspecting gauge transformation rules 
@, the first-order metric perturbation h a b is decomposed as h a b =: Ti a b + £xgab, where 7i a b and X a 
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are transformed as yH a b — xHab — 0, and yX a — xX a = C(i)o under the gauge transformation @, 
respectively.^. The gauge invariant part Ti. a b of h a b is given in the form 

H ab = ~2a 2 ^ {di 1 ) a {di 1 )b + 2a 2 v] {dri) (a {dx l ) b) + a 2 ^-2 * 7y+ £^ (d^) (dx J ') 6 , (3) 

■ (i) (i) ^ • (i) 
where D l i/j = Xfo'l = x\ = D % \ij — and D % := j l3 Dj is the covariant derivative associate with 

the metric 7^. In the cosmological perturbations [5], {<!>, W, and Xij are called the scalar-, vector- 
, and tensor-modes, respectively. We have to note that we used the existence of the Green functions 
A -1 =: (D l Di)~ 1 , (A + 2K)~ 1 , and (A + 3-fT) -1 to accomplish the above decomposition of h a b- 

As shown in Ref.[2], through the above variables X a and h a bi the second order metric perturbation 
l a b is decomposed as l a b =: C a b + c 2£xhab + \£y — £ X ) 9ab The variables C a b and Y a are the gauge 
invariant and variant parts of l a b, respectively. The vector field Y a is transformed as yY a — xY a = £? 2 ) 
+ [^(!),X] a under the gauge transformations The components of C ab are given by 



£-ab 



-2a 2 $ (dr)) a (dr]) b + 2a 2 v] (d»j) (o (da; < ) 6) + a 2 f -2 * 7y+ Xtf^) (da*),^)*. (4) 



(2) (2) ( 2 ) . (2) 

where D* z/j = xry] = X l j = f Xij = 0. As shown in Ref. |2 , by using the above variables X a and Y" a , 
we can find the gauge invariant variables for the perturbations of an arbitrary field as 

(1) Q:= (1) Q-£xQo,, (2) Q:= (2) Q-2£ x ^Q-{£ Y -£ 2 x }Q a . (5) 

As the matter contents, in this article, we consider a perfect fluid whose energy-momentum tensor is 
given by T a b = (e + p) u a u b + pS a b . We expand these fluid components e, p, and u a as 

, W 1,9 (2) , (!) l,n (2) , (1) 1 , (2) 

e = e + Ae+-A 2 e, p = p + Xp+-\ 2 P, u a = u a + A u a +-A 2 it a p. (6) 

Following the definitions ([5]), we easily obtain the corresponding gauge invariant variables for these 
perturbations of the fluid components: 

(1) (1) (1) (1) (!) W, (2) (2) (1) „ 

£ := e -fje, 7?:=P -£xp, U a :={u a ) -£ x u a , £ := e -2£ x e -{£y-i; x }e, 

(2) (2) (1) (2) (2) n) 

7? := P -2£ x P W„:=(u„) -2£ x ui-{£y-4}ii a . 

, , (1) (2) 

Through g u a Ub = g u a u b — — 1, the components of £Y a and t/ a are given by 

(1) (1) / (1) (2) (2) / (2) (2)\ 

U a = -a $ (dr)) a + a I A v +VA (dx l ) a , U a =U n {dr]) a + a I A " + V 4 I (aV) a , (7) 

(2) f /(1)\ 2 (2) / (1) (1) / . (1) M «\ 1 

W„:= a Uj - $ - I A « + Vi - W ) I -D « + V 1 - z/ H (8) 

(1) (2) 

where D l V l = D l V, = 0. 

We also expand the Einstein tensor as G b = G b + X^G a b + \\ 2 ^G b . From the decomposition of 
the first- and the second-order metric perturbation into gauge-invariant parts and gauge-variant parts, 
each order perturbation of the Einstein tensor is given by 

{1) G a b = {1) g a b [H] + £ x G a b , ^G b = Wg a » [L] + ( 2 >C [H,H] + 2£ x ^G a b + {£ Y £ 2 X ) G a b (9) 

as expected from Eqs. (0). Here, (1) £ a b [H] and {1) G a b [C] + (2) £ a b [H, H] are gauge invariant parts of the 
first- and the second- order perturbations of the Einstein tensor, respectively. On the other hand, the 
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energy momentum tensor of the perfect fluid is also expanded as T a b = T a b + A' 1 ^T a fc + iA 2 ( 2 ^T a b and 
^T a b and ^ 2 ^T a b are also given in the form 

«T/ = «T a b + £ x T a b , < 2 >T Q fc = < 2 >r a » + 2£ x ^ a b + {£ Y £ 2 X ) T b (10) 

through the definitions ([7]) of the gauge invariant variables of the fluid components. Here, b and 
( 2 ^T a b are gauge invariant part of the first- and the second-order perturbations of the energy momentum 
tensor, respectively. Then, the first- and the second-order perturbations of the Einstein equation are 
necessarily given in term of gauge invariant variables: 

^g b [H\ = 8wG^T a b , ^g a b [C} + ^g a b [H,H}^8nG ^T a b . (11) 

In the single perfect fluid case, the traceless scalar part of the spatial component of the first equation 
(i) (i) 

in Eq. (fTTj) yields v]/ = <j> due to the absence of the anisotropic stress in the first order perturbation of the 
energy momentum tensor and the other components of Eq. pip give well-known equations [5] . We show 
the expression of the second-order perturbations of the Einstein equation after imposing these first-order 
perturbations of the Einstein equations. Though we have derived all components of the second equation 
in Eq. (jTTJ) , we only show their scalar parts of it for simplicity: 

(2) (2) 3 / .1 

AnGa 2 £ = {-md n + A + 3K - m 2 ) $ -r + - A- 1 D i D j T i 3 - -T k 



2 Hd v (A + ZK)- 1 ( A 1 /)'/),! , ' - ^r fe fe ) , (12) 
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n 

(2) (2) 



8irGa?(e+p)Di V = -2ft,A $ -2WA $ +D i A~ 1 D k T k 

-3ft,A (A + 3K)- 1 (a 1 iy DjYr' - l -T k k ^ , (13) 

(2) (2) 1 

AnGa 2 V = (ft 2 , + 3Hd n -K + 2d v H + H 2 ) $ --A- 1 D > D j T i 3 

~ (ft 2 + 2Hd n ) (A + 3K)- 1 (a-'DWjT^ - \rA , (14) 
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= ^(A + si^r^A-^Ar/-^). 

where H :— d v a/a. Tq, and in Eqs. (fT2 |) - (fl~5|l are defined by 

T := +8nGa 2 (e + p) A ( v D l $ -3A $ D k $ -8 $ A $ -3 (d v $ ^ - 12 (K + H 2 ) 



(15) 



2 J 



-4 ( 9 n A $ +WA $ ) V i -2HD k <&v k +8vrGa 2 (e + p) V t V* + -D k v{ £> (fe I/O +3H 2 v k v k 



+AA $x"= -2HD k v l xL ~7;D k v l ft X / k 



+\d, £ ft x M +n £ ft x' fc -^A jffi. D k x ml +\ Dk HI dv x k]m -\ x lm (A - K) iZ, 

, /(!) (!)\ (1) (!) (!) (!) „ (!) „ « 

-167rGa 2 [ £ + V A V +12H $ A $ -4 $ ft, A $ -4ft, $ A $ 



, (!) , (1) (1) , (1)(1) (1)(1) (1) (1) (1)(1) 



-m v j A S +2^ $ ft, S -2ft,D J *£ 

ni (1) (1) ni (1) fii 1 (1) fn (1) fii 

+2D k D {i ^ l]X km +2D^ k A Xfc, +2X - A ^ --9, X ^ A Xfei +2 ft,Aj Xflfc, (16) 
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(1) (1) t 1 ) W C 1 ) 

r 2j := 167rGa 2 (e +p) A V A V -4A $ A $ -8 $ A D 3 $ 

f (i) , (i) (i) (i) / (i)\ 2 (i) (i) , „. / 

<^ 6 A $ D k $ +8 $ A $ +2 ( d v $ J + 16W $ 9„ $ +8 (2fyW + if + U 2 ) ( 

m(i) (l) (i) (l)m / (i)(V (l)^ 1 , 

5a 2 (e+p)Z> ( i vV s) -49„ $ !>(< V) +^ n D (i + 49, A $^ fc +4WD* $f' 




+327rGa 2 (e+p)A; « 

(i)(i) (1) fii rii (1) (1) a) m m 

+ 167rGa 2 (e + p) V t Vj -2 z/ fc D k D^ i>/) +2 W A A v +A f A ^fe W A W 

/ m (1) m (1) m (1) m (1) m (1) \ 

+ I -A $ # [fc ^ -A # i?' £ ^ -2 fii A v k -4dr,n { v\v k +6H 2 V k v k J ^ 

„ tW (1) t C 1 ) (1) W(l) (!) . (1) 

-4Hd n $Xij -2d 2 $ X ij -4I? fc $ D (i x,) fe +4£> fc $ A Xij -8# +4 $ A Xij 

t (1) (!) . (1) d) w ( V t 

-4D*D (< $Xj)fe +2A $ X ii +2AA $X ife 7« 

-2A 9, Xj)k -2 9,D (i Xj)fc +2 9, A Xij +-D v d n xik Hj 

(i) (1) m (1) 1 fit (1) (1) fn (1) (i) 

+d, xl d n x/ +2D p Xfc], X/ ~2 A X» A X* fe - *' m A A Xw +2 x' m AA.i Xj) m 

(1) fii / 3 m (1) 3 m (1) l a) (1) m (1) \ 

- x lm An A x« + I x/i ft, x H + 2 A xiL D k x rnl ~2 Dh xii D l x mk +k xiLx lm J m, 

and := 7 J r^.. These equations ([T^ -lflS" ]) coincide with the equations derived in Refs.[3] except for 
the definition of the source terms To, Ti, and . Further, as shown in Refs.[3], the equations (|12[) and 

. . . (2) 

(| 15(1 are reduced to the single equation for $ . We also derived the similar equations in the case where 
the matter content of the universe is a single scalar field[4]. 

In summary, we have extended our formulation without ignoring the first-order vector- and tensor- 
modes. As the result, these equations imply that any types of mode-coupling arise due to the second-order 
effects of the Einstein equations, in principle. In some inflationary scenario, the tensor mode are also 
generated by the quantum fluctuations. This extension will be useful to clarify the evolution of the second 
order perturbation in the existence of the first-order tensor-mode. Further, to apply this formulation to 
clarify the non-linear effects in CMB physics [BJ, we have to extend our formulation to multi- field system 
and to the Einstein-Boltzmann system. These extensions will be one of our future works. 

References 

[1] C.L. Bennett et al., Astrophys. J. Suppl. Ser. 148, (2003), 1. 

[2] K. Nakamura, Prog. Theor. Phys. 110 (2003), 723; K. Nakamura, Prog. Theor. Phys. 113 (2005), 
481. 

[3] K. Nakamura, Phys. Rev. D 74 (2006), 101301(R); K. Nakamura, Prog. Theor. Phys. 117 (2007), 
17. 

[4] K. Nakamura, in preparation. 

[5] J. M. Bardeen, Phys. Rev. D 22 (1980), 1882; H. Kodama and M. Sasaki, Prog. Theor. Phys. Suppl. 
No.78, 1 (1984); V. F. Mukhanov, H. A. Feldman and R. H. Brandenberger, Phys. Rep. 215, 203 
(1992). 

[6] N. Bartolo, S. Matarrese and A. Riotto, JCAP 0401, 003 (2004); N. Bartolo, S. Matarrese and 
A. Riotto, Phys. Rev. Lett. 93 (2004), 231301; N. Bartolo, E. Komatsu , S. Matarrese and A. Ri otto, 
Phys. Rept. 402, 103 (2004); N. Bartolo, S. Matarrese, and A. Riotto, |arXiv:astro-ph/0512481| 



4 



